On the Deconfinement Phase Transition in Hot Gauge Theories with 

Dynamical Matter Fields 

O. Borisenko, M. Faber 

Institut fiir Kernphysik, Technische Universitdt Wien, 
Wiedner Hauptstr. 8-10, A-IO4O Vienna, Austria 

G. Zinovjev 

Institute for Theoretical Physics, National Academy of Sciences of Ukraine, Kiev 

252143, Ukraine 



Abstract 

The phase structure of hot gauge theories with dynamical matter fields 
is reexamined in the canonical ensemble with respect to triality. Since this 
ensemble implies a projection to the zero triality sector of the theory we in- 
troduce a proper quantity which is able to reveal a critical behaviour of the 
theory with fundamental quarks. We discuss the properties of both the chro- 
moelectric and chromomagnetic sectors of the theory and show while electric 
charges carrying a unit of Z{Nc) charge are screened at high temperatures 
by dynamical matter loops, this is not the case for the Z{Nc) magnetic fiux. 
An order parameter is constructed to probe the realization of local discrete 
Z{Nc) symmetry in the magnetic sector. We argue it can be used to detect 
a deconfinement phase being defined in terms of the screening mechanism as 
a phase of unscreened Z{Nc) flux. It may be detectable at long range via the 
Aharonov-Bohm effect. We discuss the possible phase structure of QCD in 
this approach. 
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1 Introduction 



A new round of considerable interest in discrete global and local gauge symmetries 
started when it was conjectured that the discrete center Z{Nc) of an underlying 
gauge group G can be of crucial importance for quark confinement |jl], 0]. Another 
mile-stone appeared with understanding the deeply rooted relation between the 
spontaneous breaking of the global Z{N^ symmetry and the deconfinement phase 
transition in pure gauge models on the lattice and in continuum |0. This paper 
continues the previous studies of hot gauge theories with dynamical matter fields 
and emphasizes the importance of the magnetic sector in the investigation of their 
phase structure reliably. We propose a new order parameter to test the realization 
of the discrete Z{Nc) symmetry at finite temperature and to measure screening 
effects in different regions of temperature and couplings. The result seems equally 
applicable to any Z{Nc) and SU{Nc) gauge models coupled to Higgs and/or fermion 
fields carrying Z{Nc) charge, i.e. which are nontrivial on the center group. 

A brief overview of commonly recognized results of finite temperature QCD 
analysis occurs quite pertinent for an introduction to the subject. So then, gluon 
fields are strictly periodic in time, while quarks are antiperiodic with a period given 
by the inverse temperature. Pure gauge theory has an exact Z{Nc) global symmetry. 
The gauge invariant operator, the Polyakov loop (PL), = -^TrYl^^^Uo^Xjt) 
transforms under Z{Nc) global transformations as 

2'jii 

ZL^, Z = exp[— n], n = 0, - 1. (1) 

The PL can be used as an order parameter to test Z{Nc) symmetry in pure gauge 
theory. The expectation value of the PL is interpreted as the free energy of a probe 
quark Fq immersed in a pure gluonic bath 

<L^>=exp(-^F,). (2) 

Unbroken Z{Nc) symmetry implies < L >= and Fq = oo. When the global 
Z{Nc) symmetry is spontaneously broken, < L > develops a non-zero value giving 
a finite magnitude to Fq, i.e. it costs a finite energy only to create a single quark 
in the gluonic bath. However, in it was discussed that in a system with a 
finite ultraviolet (UV) cutoff the free energy should not diverge. Usually, Monte 
Carlo (MC) and analytical calculations are also performed with periodic boundary 
conditions (PBC) in space directions. It has been shown, however, that the Gauss' 
law and this PBC in space are inconsistent unless the sum of quark and gluon colour 
charges vanishes, what looks controversial since they have different triality. In fact. 



2 



it was concluded that for space PBC the expectation value of the PL is not the free 
energy of heavy quark. 

The second trouble appears if we realize that in the spontaneously broken phase 
< L > may pick up A^^^ different values corresponding to A*",, equivalent minima 
of the free energy. Thus, L can be negative or even complex and Eq.(|D tells us 
that the free energy could be some complex number. This gives rise to doubts 
that Eq.(^ has the proper physical interpretation. When dynamical quarks are 
included the picture becomes more comphcated and new troubles appear. The 
fermion determinant generates loops winding around the lattice a number of times 
which is not a multiple of Nc- Such loops present a propagation of single quarks 
and transform non trivially under Z{Nc). This means that dynamical quarks break 
Z{Nc) symmetry explicitly and screen sources of heavy quarks at any temperature. 
The expectation value of the PL prefers the phase with argL = which provides the 
minimum of the free energy. Other Z{Nc) phases with argL = k = 1, . . . , N^ — l 
become metastable. They possess, however, such unphysical properties as a complex 
free energy or entropy Recently, it has been discovered that chiral symmetry is 
not restored in Z{Nc) phases 0, what led to another reexamination of degenerate 
Z{Nc) phases and interfaces between them in pure gauge theory and one of the 
conclusions is that all Z{Nc) phases correspond to the same physical state and the 
interfaces are unphysical. The operator, which we introduce here may test directly 
the Z{Nc) phases and hence can be used to clarify this important problem in the 
description of finite temperature QCD. 

The main motivation of our study comes, however, from the following. It became 
almost a dogma to consider that the deconfinement phase transition is related to 
the appearance of nonzero triality above the critical point. A number of 'order 
parameters' were constructed aiming to detect a phase transition and nonzero triality 
states 1^. To our knowledge all of them failed to display a phase transition and 
to test any nontrivial triality state. The most conclusive results here have been 
obtained in Z{Nc) gauge theories with Higgs fields p. Since it had been known that 
at zero temperature this system has two phases, a confining/ screening phase and a 
deconfining phase [1^ separated by a critical line and since such a critical line was not 
found at finite temperature it was concluded that presumably at finite temperature 
the critical behaviour is not present at all. Moreover, it has been claimed that there 
is no smooth relation of finite and zero temperature theory because no phase has 
been found at arbitrary small (but finite) temperature which could correspond to a 
deconfining phase at zero temperature in the same interval of bare coupling values. 
By now it became almost usual to define deconfinement as a 'crossover' but not a 
genuine phase transition in full theory and to study the 'phase transition' related 
to chiral symmetry restoration at light quark masses |TT|. The question how this 
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phase transition is related to the deconfining crossover at larger quark masses has 
not been answered positively yet. 

Two points are to be discussed here. First of all, the deconfined phase was 
identified with the phase where nonzero triality should exist and order parameters 
have been constructed to detect such 'free triality states'. The correctness of such 
a conjecture is not obvious a priori. In the context of this model, we demonstrate 
here that there might be a phase transition which, however, is not related to the 
liberation of triality but rather to different screening mechanisms of triality and 
introduce an operator which tests these different mechanisms. Second, the previous 
emphasizes have been put on the possible realization of the Z{Nc) global symmetry 
in the electric sector of the theory as explained above. The magnetic sector has 
been untouched from this point of view. We intend to overcome this deficiency in 
the present paper putting the problem on a rigorous basis and suggesting to study 
local Z{Nc) symmetry in this sector. 

Certainly, we need an appropriate framework to study the problem. The grand 
canonical ensemble (GCE) with respect to triality is obviously not a proper tool 
since the Z{Nc) symmetry is explicitly broken in this ensemble. In our previous 



papers we introduced an ensemble canonical with respect to triality [12| (see also 



1^) to analyze the QCD phase structure [|T^, |T5|. Since the Z{Nc) symmetry is not 
explicitly broken in the canonical ensemble (CE) it is more suitable to study the 
QCD phase structure, especially the chromoelectric sector, therefore we shall use 
this formulation in the present paper as well. Of course, in general we expect in the 
thermodynamic limit the same phase structure both in the GCE and in the CE. 
The paper is organized as follows. 

In Section 2 we introduce the Z{Nc) invariant quantities. At and Ag, which can 
probe a unit of Z{Nc) charge or fiux and discuss its properties. These quantities 
give an opportunity to treat the chromoelectric and chromomagnetic sectors of the 
theory on the same footing applying them for the investigation of the ground state. 
In Section 3 we discuss in some details the general phase structure of gauge theories 
at zero and finite temperatures and outline a realization of local Z{Nc) symmetry 
in the chromomagnetic sector of the theory different from the usually discussed 
chromoelectric sector. We define the deconfinement phase by analogy with the 
zero temperature deconfinement phase, phase where the kinetic screening 

due to the gluon-gluon interaction is stronger than the dynamic screening due to 
dynamical matter loops and propose an order parameter to test the corresponding 
phase transition. Further, we explain how to investigate the problem of domain 
walls in QCD and suggest a possible resolution. To check our predictions we study 
in the Section 4 a model of Z{Nc) gauge fields coupled to discrete Higgs fields. 
We calculate the order parameters At and Ag in different sectors and for various 
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regimes of coupling constants and temperature. Finally we summarize our results 
and outline perspectives for the application of our results for Z{Nc) and SU{Nc) 
models with dynamical quarks. 



2 An order parameter for Z{Nc) charges and do- 
main walls 

In this section we construct an order parameter at finite temperature which may 
probe a nontrivial charge and/or flux in the Debye screened phase and can be rep- 
resented by a Z{Nc) invariant operator in a finite volume. Such a quantity was 
introduced some years ago in the context of gauge theories with nontrivial discrete 



abelian center at zero temperature |16]. We will discuss that this quantity is able 



to detect a nontrivial phase of some state in a region enclosed by a surface S 

with the limits Tj,C oo taken. W{C) is a space-time Wilson loop nontrivial 
on the Z{Nc) subgroup. S may be thought as a set of plaquettes in a given time- 
slice whose dual plaquettes S* form a closed two-dimensional surface for this time 
slice. Further, E has to be chosen in such a way that it encloses at the considered 
slice a time-like line of the loop C. As will be discussed below, the operator -F(E) 
frustrates the plaquettes of E by the so-called "singular gauge transformations" with 
a non-trivial element of Z{N). 

The action of this order parameter is based on the Aharonov-Bohm effect: despite 
the absence of electric fields a singular potential influences particles at arbitrary long 
distances giving them a nontrivial phase during winding around the solenoid. In our 
case a Z{Nc) charge may be viewed as a kind of such a singular solenoid placed at 
some space position. To get this configuration one should take the Wilson loop 
W{C) and consider the limit C —>■ oo. Since the Wilson loop is a source of Z{Nc) 
charge one can probe it via the "device" -F(E). A nontrivial Z{N^ charges may be 
detected only in the case when the fermionic screening is suppressed relatively to 
the Debye screening 



lim A(E, C) = exp 



(4) 



where K{J], C) is the linking number of the surface E and the loop C. If the triality 
charge is totally screened by fermions one gets lims,c^oo ^(S, C) = 1. 
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Let us discuss the implementation of this order parameter for gauge theories 
at finite temperature. The following discussion is equally applicable for any gauge 
theory if it admits an introduction of the PL which transforms nontrivially under 
the discrete center of the gauge group, i.e. SU{Nc) or Z{Nc). We consider the 
lattice gauge model with a standard Wilson action given by 

Sw = XT.TrUp, (5) 

p 

where Up is a product of gauge link matrices around a plaquette p. Then, since the 
fundamental PL is a source of Z{Nc) charge at finite temperature one can insert the 
correlation function of PLs into @ instead of a Wilson loop to get 

The operator F which frustrates the plaquettes of S is defined by 



F(E) = exp 



Xj2{Z-l)TrUp 

pes 



(7) 



where Z is a non-trivial element of Z{Nc). One can then define the "frustrated 
action" in the following way 

^ = A ^ TrUp + ZTrUp . (8) 
p^E pes 

A schematic diagram of the surface S in space is shown in Fig.l for the case of a 
(2 + l)-dimensional theory. The two Polyakov loops Lq and Lr are closed in time 
direction. The frustrated plaquettes S are drawn, they enclose the loop Lr. The 
links dual to the plaquettes of S form a closed line S* surrounding Lr. They are 
shown by a dotted line. The frustration of the plaquettes of S can also be achieved 
by multiplying the time-like links enclosed by S with Z. Bold lines represent these 
links. 

The canonical partition function reads 

Z = ^lljX{dUiX{ d¥jmi e-^--^'^'^) . (9) 

^ k=l I x,i 

Sq is the standard quark action, where we have to substitute Uq — > exp[^^]f/o for 
one time slice in the CE |T2[. Then, J2k projects to zero triality states. We define a 
new partition function in the CE with the action Sp (p|) as 



1 iVc ^ 

W E / UdUiUdKdK e-^--^^(^') . (10) 
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In what follows we refer to the expectation values < ■ ■ ■ >f calculated in the ensem- 



ble defined by (10) as those obtained in the F-ensemble to distinguish them from 
the corresponding values < ■ ■ ■ >o in the standard ensemble. It is easily seen that 
with such a definition At in reads 



< LqLr >o 




Figure 1: Correlation function of Polyakov loops in the F-ensemble for the case of 
a (2+l)-dimensional theory. The two Polyakov loops Lq and extend in time 
direction. The frustrated plaquettes are explicitly depicted. Their dual links form a 
closed (dotted) line around Lr. Frustration can be achieved by multiplying of the 
bold links with a nontrivial element Z G Z{Nc). 

To probe a Z{Nc) charge in space one has to take the limit R ^ oo where 
correlations of the ratio (^) reduce to one PL. Thus, one can probe a unit of Z{Nc) 
charge enclosed by the surface S. 

A similar F-ensemble can be constructed for the spatial Wilson loop. In this 
case the order parameter A is essentially the same as in zero temperature theory, 
i.e. it is defined in formula (§) where we should insert a pure space-like Wilson loop 
instead of a time-like one. In what follows we use the notation Ag to distinguish 
it from At introduced above. The meaning of the operator in this case is however 
different as we shall discuss in the next section. 
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Let us investigate now some simple properties of At and Ag. In a pure gauge 
system the surface S may be completely gauged away by a change of variables 
Uo{x) Z*Uq{x) on the time-like links in the three dimensional volume enclosed 
by the surface S. It follows that A = Z ^ 1 for the pure gauge system since 
the fundamental PL in the origin changes sign, precisely as in the theory at zero 



temperature [T^. When dynamical fermion (Higgs) fields are added, the surface E 
cannot be gauged away since the corresponding time-like fields Uq will be flipped in 
the fermionic (Higgs) part of the action, so that 

S = XY,TrUp + a ^ Qi{Ui) + a QiiZ*Ui) , (12) 

where Qi means the density of either the fermion or the Higgs action and f2(S) is 
a volume enclosed by the surface S. X)/ is the sum over links. The same obviously 
holds for As with the corresponding substitution of space-like links instead of time- 
like ones. 

An important interpretation for the influence of -F(S) comes from the following 
observation. -F(S) tries to implement a phase change at the spatial surface S. If the 
global Z{Nc) symmetry is spontaneously broken -F(S) produces a stable interface 
between the volume enclosed by S and the surrounding vacuum. The different phases 
can be detected by the PL and the order parameter At{J],C). If the symmetry is 
unbroken F(S) cannot produce a stable interface and the value of the PL is not 
influenced by F{T,) which is far apart. 

The most essential property of At and As is, however, that they measure screen- 
ing effects of dynamical fermion (Higgs) loops (dynamical screening) together with 
screening effects of the kinetic energy due to the pure gluonic interaction and show 
which screening is stronger at given conditions. The 'Aharonov-Bohm effect' works 
only when screening due to fermion (Higgs) loops is weaker so that it is possible to 
detect a charge (flux) at spatial infinity. 

The simplest example where the quantity A can be calculated exactly is the 
Gaussian model for a scalar field coupled to an external field. We want now to 
consider this example which may be quite suggestive since it illustrates how the 
operator works; in this case it does not correspond to any discrete charge but rather 
shows the competition between two kinds of screening and measures which one is 
stronger. Let us define the partition function of the massive Gaussian model in an 
external field h^j. as 



Zo=[ ]]_da^exp (-^-'^{ax - a^+nf - m^'^fTl + 'Yh^orA . (13) 

X \ ^ x,n X X ) 

In this case S consist of links dual to a (D — l)-dimensional surface S*. The operator 
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F(E) is introduced on the stack of links E and has the form 

= exp 





= exp 


~2 ^ ax(Tx+n 






/6E 



(14) 

In what follows we consider a D-dimensional periodic lattice with period L and 
number of sites V — LP . The region Q enclosed by E we choose to be a — 1 torus 
winding around the lattice. Performing now the substitution cr — > — cr in the interior 
we can write down the partition function in the F-ensemble as 



/oo ^ 
exp 
-oo ™ 



a;^Q(S) a;GQ(S) 



(15) 



Figure 2: Correlation function of spins in the F-cnscmble for the case of a 2D theory. 
The frustrated links E arc explicitly depicted. Their dual links form two closed lines 
around and are chosen to wind around the whole lattice. 

After this transformation we have 

^^_ <a,an >f ^ ^^^^ 
< ctqCTr >0 

where < (jQap >o and < aoa^i >p are the correlation functions caknilatcd in the 
usual and the F-ensemble, correspondingly. A diagram of the correlation function 
in the F-ensemble is plotted in Fig. 2. On a finite lattice < aoaR >q is known to be 

< CToaR >o 



/L-l 



2V 



1 .2TiRk^ 
+ 4^ (^E^fe ^fcexp[z— ^] 



yk=o / 



(17) 
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where 

hk = -j= 2^ Kexp[t——\ (18) 

V V/ fe=0 ^ 
and the propagator Mj7^ in the momentum space has the well known form 

M^' = {D + m^-Y.cos^)-\ (19) 



For < (jQaji >F we should change signs for in the Fourier transform ([T8| ) at 
corresponding sites. Let us consider now the simple case of a constant external 
field. We have from (|T8|) hj. = \/VSk,oh. In the thermodynamic limit V ^ oo we 
get the following expression 



G{R) + hy 2 (m^' - 2f{R)) (M q-^ - 2/(0)) 
G{R) + hy2{M^^y 



A=- ' , j:^:!': — (20) 



where we used 

GiR)=l (^1 ^-—T^^T^r^——. (21) 




m) = E / ? TT^^ 7 (22) 

JO 27r 1 + m"^ — COS0 

and Mq = im? . Lp in the last formula denotes the number of sites with negative 
sign of h. In any dimension G{R) goes to zero exponentially as R increases. In the 
limit R and Lp ^ oo one can get the final expression for A 

^=^^. (23) 



V2 



m 



One can see from the last formula that for a constant external field A is always 
positive. It is interesting to note that the limit m exists even for the one- and 
two-dimensional cases and gives A = 0. The limit m oo corresponds to the value 
A = +1. This is the only mass when A = 1 since the scalar field a does not carry a 
unit of any discrete charge. As we discussed earlier in Eq. (^, for fields carrying a 
discrete charge we would find only discrete values of A. 

The interpretation of these results is rather transparent. There is a competition 
in the correlation function between the term G{R) and the second term of (|T7|). G{R) 
comes from the kinetic energy of the system, decays exponentially and describes the 
corresponding Debye screening. The second term comes from the external field and 
describes "screening" due to this field. We think this example shows that A is able 
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to determine which screening mechanism is stronger. In fact, it is not difficult to find 
a nonconstant external field such that in the limit R ^ oo the quantity G{R) goes 
to zero slower than the second term of the correlation function and hence A = —1. 
A special example of this type we will mention later. 

3 Screening at zero and finite temperatures and 
the deconfinement phase transition 

We are now ready to proceed to study the phase structure of gauge models with a 
nontrivial center. We would like to explore the above mentioned properties of At 
and As to reveal some features of these models. Let us start with a discussion of 
pure gauge theory at zero temperature. First, we are interested in Z[Nc) gauge 
models. These models are known to exhibit confinement at strong coupling and 
deconfinement behaviour at weak coupling. This means that the Wilson loop obeys 
an area law in the first case and a perimeter law in the second. The interpretation 
of such a behaviour in terms of the quantity A has been given in [T^ and amounts 
essentially to: A = —1 and the absence of dynamical screening of a Z{Nc) charge at 
any coupling. The Wilson loop in the F-ensemble always changes its sign, though in 
different ways, in strong and weak coupling regimes. At weak coupling the ground 
state of the system consists of gauge spins flipped inside a volume bounded by 
S relatively to gauge spins outside. A perimeter law for the Wilson loop arises 
in this phase due to the kinetic energy of gluons as can be seen most easily in the 
Hamiltonian formulation (see [0 for a comprehensive discussion of this topic). If we 
add Higgs flelds to this system dynamical screening appears, causing a competition 
with the screening coming from the gluon kinetic energy. Before we proceed further 
we need to give more precise deflnitions of dynamical screening and of screening 
coming from the gluon kinetic energy. 

Let W{C) be a Wilson loop in Z{Nc) lattice gauge theory either at zero or at 
finite temperature. In the weak coupling phase of pure gauge theory one may write 

<W{C)>^exp{-^giP), (24) 

where P means the perimeter of the loop C. We refer to nonzero 'jgi in the decon- 
finement phase as coming from the kinetic energy of the pure gluonic interaction 
and attribute the corresponding kinetic screening to this interaction. Let us imagine 
now this system coupled to dynamical matter fields, either to Higgs or to fermion 
ones. Among other contributions to the Wilson loop the dynamical fields generate 
a perimeter decay of the loop at any coupling and/or temperature. In the C ^ oo 
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limit one always finds the following behaviour of the loop in the strong coupling 
limit, i.e. A = in (^) 

<l^(C)>~exp(-7,,„P). (25) 

We refer to '^dyn as coming from the dynamical screening due to the dynamical Higgs 
or fermion fields. 

Let us suppose now that the system is in a region of coupling constants where 
the pure gauge interaction produces only area law decay for the Wilson loop, i.e. 
A << 1. One has the following formal expansion 

< WiC) >~ exp(-7rf,„P) + K2exp{-agiS) + ... , (26) 

where agi is a string tension of the pure gauge theory and S is the minimal area 
enclosed by the loop C. Dots mean all the (perimeter and area) terms which vanish 
faster in the thermodynamic and C — *■ cxo limits than the corresponding written 
terms. From the properties of A described in the previous section it is clear that in 
the F-ensemble the term with area law decay should change its sign. The expansion 
might look like 

< W{C) >~ iriexp(-7,,„P) - K^expi-agiS) + ... . (27) 

Substituting these expansions into (^ and taking the limit C — > cxd we find A = 1 
since the area term coming from the pure gluon interaction vanishes clearly faster 

Thus, A takes a trivial value and we interpret this behaviour as reflecting the follow- 
ing fact: in the confinement region the dynamical screening is the only one screening 
present in the system. Actually, in this region there is not even a real competition. 
However, in the pure gauge system there is a critical point which of course is not a 
critical point of a full system. Above that one has instead of ( p6|) 



<W{C)>r^Ksexp{-jay^P)+K^exp{-jgiP) + ... . (29) 

In the F-ensemble one has to change the sign of K4. There is a competition between 
the kinetic screening and the dynamical screening. One gets 

For the second possibility the kinetic screening is stronger and a Z{Nc) charge can 
be detected. We interpret this as a permanent feature of the deconfinement phase. 
Simple and maybe well known conclusions are: 
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1) deconfinement is the phase at a small gauge coupling constant; 

2) deconfinement may take place even in the presence of dynamical screening; 

3) A acts equally good for time and space-like Wilson loops since there is no 
actual difference between them in zero temperature theory. It distinguishes between 
confined/screened phase and deconfined phase. 

The first observation is rather important especially in the context of finite tem- 
perature theory. Usually, the deconfinement phase is associated with the high tem- 
perature region identifying the latter with the region of small couplings. In fact, at 
finite temperature the original formulation involves two free and completely inde- 
pendent parameters: the temperature and the bare coupling. The region of small 
bare coupling needs a priori not to be the phase of high temperature. In fact, rig- 



orous results on the deconfinement transition in pure gauge models [ITB] state that 
as soon as the temperature is turned on there exists such a coupling below that the 
system is in the deconfined phase. 

It is also interesting to mention that the equation (|30| ) predicts an exact equation 
for the critical line in the theory, namely 

Gdynia) = IglW , (31) 

where A and a are the gauge and Higgs (fermion) couplings, respectively. 

The main question arising from these observations which we would like to put 
here concerns special features of the deconfined phase at finite temperature relatively 
to the one at zero T? We argue that in fact there are no special features, breakdown 
of global and/or local Z{N) symmetry at high temperature may be misleading in this 
respect if one deals with the full theory. Triality is always a conserved quantity and to 
reveal the critical behaviour one has to look at the screening mechanisms of triality in 
different regions of coupling. We therefore suggest to define the deconfinement phase 
in the theory with dynamical matter fields as the weak coupling phase where the 
kinetic screening due to the gluon interaction is stronger than dynamical screening. 



A similar proposition was discussed also in [|T3|. Technically however the situation 
at finite temperature may differ from the picture discussed above. 

Let us discuss now in detail what we expect for finite temperature theory. We 
continue considering Z[Nc) gauge models in the Euclidean path integral formula- 
tion. There is no symmetry anymore between electric and magnetic sectors at finite 
temperature, in particular, one can find a drastic difference in the behaviour of the 
corresponding Wilson loops, i.e. the spatial Wilson loop and the thermal (Polyakov) 
loop. The spatial Wilson loop behaves essentially as in the zero temperature case. 
Thus, the first observation which can be easily deduced from this fact is that the 
quantity Ag should be as good an order parameter at finite temperature as it is at 
zero temperature because all formulae above are valid in this case if we substitute 



13 



W{C) by the spatial Wilson loop and A hj As. A nontrivial value of As should 
imply that a unit of Z{Nc) flux is not screened dynamically and can be detected at 
long-range. Therefore, if there is a critical behaviour at zero temperature it should 
also be present at finite temperature. We would like to stress that this observation 
does not depend on the ensemble (CE or GCE) used to study the finite temperature 
system (see also next section). 

On the other hand, the behaviour of order parameter At could be qualitatively 
described as follows. In the strong coupling phase the correlation function of PLs 
calculated in the pure gluonic sector goes down exponentially with increasing dis- 
tance. The fermionic sector, however generates terms which screen heavy quarks 
and lead thus to a constant value of the correlator even at spatial infinity. This 
implies At = 1. In the weak coupling region the pure gluonic sector also gives fi- 
nite values for the correlation at spatial infinity leading in such a manner to the 
competition with dynamic screening. In such a situation the direct use of At as 
an indicator of a phase transition is not possible because both contributions stay 
finite in the i? — > cxd limit. One may argue, however that At = 1. To understand 
qualitatively what one should expect in such a situation one can study some simpler 
models with similar properties. Such studies were done by us in |0 on the example 
on the 3D Ising model where we were able to define and to calculate an analogue 
of the quantity A for spin systems in the presence of an external magnetic field. 
It was done both analytically and in MC simulations. Applying that consideration 
to the Z{2) gauge model we conjecture that in the weak coupling phase the value 
At = 1 has the following meaning: The operator F introduces a stable interface 
in the pure gauge system. In the weak coupling region the main state is the one 
where all the PLs are flipped inside the volume enclosed by the surface E relatively 
to the loops outside of the volume. At = —1, and this means that the interfaces 
are stable. Dynamical matter loops try to destroy such interfaces. At sufficiently 
weak coupling the main state has all loops aligned in one of the Z{Nc) direction (in 
the canonical ensemble). The state with the loops flipped inside of S, which could 
potentially lead to a nontrivial value appears to be suppressed by a volume term 
whereas the contribution of the state with aligned loops is suppressed with a surface 
term. There is a necessary mathematical condition for At to take a nontrivial value: 
the expansion for large masses or small Higgs coupling has to converge in the electric 
sector and in the thermodynamic limit. The consideration in |^ shows that such an 
expansion is presumably not convergent at any finite temperature. It means in turn 
that it is impossible to expand the correlation of the PLs similarly to the Wilson 
loop at zero temperatures, e.g. (^) and (^). Indeed, if one wants to make such an 
expansion one would find at the first nontrivial order in small Higgs or fermion and 
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large inverse gauge couplings the equations for At (see next section for details) 



where A is gauge coupling, V is the volume and the constant C is either proportional 
to the Higgs or to the fermion coupling. A similar expansion one can find for the 
correlator in the F-cnsemble. One sees from here that the thermodynamic limit 
does not exist already at the first nontrivial order, even for the correlation of the 
PLs. One concludes that in this situation At cannot determine which screening is 
stronger. However, it gives an interesting information on the domain wall structure 
of the high temperature phase. The value At — 1 corresponds as explained above to 
an instability of the interfaces of pure gauge theory. Wc argue thus that the order 
parameter At can be used to resolve a long standing problem of high temperature 
QCD, namely the problem of the stability of domain walls. In the CE all Z[N^ 
directions are degenerate but in the thermodynamic limit, as the behaviour of At 
shows there is no stable interfaces between regions with different orientations of the 
PL. In the next section we mention however a possibility where the interfaces could 
be stable and At = —1. 

Let us shortly summarize. In the confined phase a Z{Nc) charge does not exist 
in the sense that quarks are bound in triality zero states, i.e. mesons and baryons. 
Whatever small it is, dynamical screening is finite at any long range. Thus, a 'free' 
Z{Nc) electric charge cannot be found and A = 1. In the weak coupling phase, 
quarks may become free and kinetic screening appears to play an important role. 
The correlation function is finite because of this screening which is defined as coming 
from the pure gluonic part. It is however rather impossible to distinguish between 
this and dynamical screening because both of them lead to a finite value of the 
correlation of the PLs at spatial infinity. The attempt to separate the different 
types of screening leads to volume divergences. This fact does not imply that the 
critical behaviour is lost, it concerns only purely time- like loops. One can use the 
quantity Ag to measure screening in the magnetic sector, i.e. the presence of a 
unit of the Z{Nc) flux. Such a treatment should reveal the critical behaviour if it 
is present at zero temperature. Indeed, the screening mechanism in the magnetic 
sector is not much affected by dynamical matter loops winding around the lattice 
in time direction and should be essentially the same (in the thermodynamic limit) 
as at zero temperature. 

In the next section we show how these ideas work in the theory of Z{2) gauge 
spins coupled to Higgs flelds. 




(32) 



15 



4 Application: A in Z(2) gauge theory 



As an application of the proposed approach we would like to examine here the 
model of Z{2) gauge spins coupled to discrete Higgs fields at finite temperature. 
The canonical partition function of the Z{2) gauge model is given by the path 
integral 

^ = ^ E E E . (33) 

<T=±1 Si=±l Za:=±l 

The sum over a is the sum over two n-ality sectors. We have denoted 

= E Wo + E^n^pn , (34) 

PO Pn 

where Up^ iUp^) is a product of Z{2) gauge link variables s^{x) around a space(time) 
plaquette. The Higgs action at finite temperature we write down as 

aso{x)zx+o , (35) 

x,n X 

where both gauge s^{x) and Higgs fields obey a periodicity condition. The phase 
diagram of this theory at zero temperature was well established long ago [IC, 1^ and 



was discussed in the previous section. The phase diagram at finite temperature is 
rather unclear. It has been recently claimed that pure Z{Nc) gauge model may well 
possess 4 distinct phases in terms of bare coupling depending on actual values of Aq, 
A„ and Nt [^. Ref. claims that the coupling of gauge spins to the Higgs field 
leads to a trivial phase structure, the system is in a confining/screening phase at all 
couplings and temperatures. The U{1) abelian Higgs model at nonzero temperature 



was studied earlier in [|21[. It has been argued that there is no phase transition in 
temperature, neither increasing the temperature from the confining/screening phase 
nor from the deconfined phase. However, there should be a phase transition at finite 
temperature in the coupling constant from the confining/Higgs phase to the decon- 
fined phase. It was speculated that this phase boundary may disappear above some 
finite temperature but the issue remains unsettled. We expect a similar scenario to 
be realized in the Z{2) Higgs model, i.e. the deconfinement phase transition is a 
transition in the coupling constant either at zero or at some finite temperature. 
At zero temperature the order parameter A was calculated for the case of Z[2) 



gauge spins coupled to Higgs fields in Ref.|[T^. It was shown that A indeed is 
an order parameter changing abruptly from A = 1 (confining/dynamical screening 
phase) to A = —1 (deconfining phase). 
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Let us study first At in the finite temperature Higgs model. Shifting the surface 
S via Z{2) singular gauge transformations to the Higgs part of the action we get 



= (36) 

where the PL is given by 

Nt 

Lo = Y[soiO,t). (37) 
t=i 

The time-like part of the Higgs action reads 

= '^ho{x)Zx(JSo{x)Z:,:+o , (38) 



where 



and fl is as in (0). We would like now to analyze two regions of the gauge cou- 
plings Ao and A„, the weak and the strong coupling regions, at arbitrary Nf. The 
corresponding expansions are known to be convergent, precisely as in the zero tem- 
perature theory [1^, . 

I. At in the strong coupling region, Aq and A„ << 1. 

We want to get an expression for the partition and the correlation functions 
which could be used in both ensembles and both for At and Ag. For that purpose 
we introduce a space-time dependence in the Higgs couplings and use the general 
notation h^. Following the standard scheme of the strong coupling expansion we 
get (up to the third order), e.g. for the partition function the expression 

Z = C{\) Y\_cosh.h^{x)[l + tanh A^ ]~[ tanh h^{x) + 
/ p i&p 

tanh A^ tanh A^/ tanh /i^(x) + ^ tanh A^ tanh A^/ ]^ tanh/i^(x)] (40) 
PT^p' i&p^p' pi^p' ^^p^p' 

and we introduced the obvious notation A^. X]p means sum both over time-like and 
space-like plaquettes, means that one has to omit plaquettes which have a link 
in common and means sum over plaquettes which have one link in common, 
njgp p/ means that one has to omit that link from the product which is common for 
both plaquettes. At last, 

C(A) = (cosh Ao)^™ (cosh A„)^^", 



17 



where Np^, Np^ is the number of time(space)-hke plaquettes, correspondingly. The 
general expansion for the correlation of the PL is 

< LoLr >= ic(A) E LoLrc^'^II + 

^ si=±lz^=±l 

J2 tanh Xf, Sp+J2 ^anh tanh A^^ SpSp' + 0(A^)]. (41) 

Now it is straightforward to calculate correlations of the PL in both the standard and 
the F-ensemble. Let us suppose that the frustrated links are in the time slice t = 1. 
Since in this case space couplings hn{x) are not affected by a singular transformation 
one can omit the space-time dependence. It gives up to second order 

Nt 

< LqLr >= tanh /io(0,t) tanh ^o(-R, if:) [1 + 2-DA^t tanh Aotanh^ /i„(l-tanh^ ho)], 
t=i 

(42) 

where D is the space dimension. Since the linking number of the PL in the origin 
and the surface S is 1, it gives the result 

A = 1 + 0{X'). 

In fact, it is obvious that the expression in the square brackets of (^) is an even 
function of hoit) up to the order (tanh Aq)''"^, where Ls is the radius of the region 
Q. Only on the boundary the corresponding plaquettes will change signs. Thus, 

A = 1- A^sCl(tanhAo)^^ (43) 

where iVs is the number of frustrated plaquettes and Ci is independent of E. It 
leads in the S ^ oo limit to the expected result At = I. 
II. At in the weak coupling region Aq , A„ >> 1. 

We fix a static gauge where all gauge spins in time direction So{x) are set to 1 
except for one time-slice t = 1 which includes the frustrated plaquettes on S. The 
Debye screening comes from the pure gluonic action and produces the formula for 
the correlation of PLs 

In < LqL/j > oc exp[-M^i?], (44) 

where Md is the Debye mass. Though it may seem at the first glance that the 
fermionic contribution to the correlation function is suppressed as /iq^* and could 
be neglected relatively to the Debye screening, this is misleading. The crucial point 
is that in order to skip this contribution one should isolate the contribution from 
the fermionic sector. To do that, one needs to make an expansion of the correlation 
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function in small Hq. But such an expansion is known to be divergent |^. Thus, 
there is no direct way to separate the Debye screening from the fermionic screening 
in the electric sector. In this situation to calculate At we should construct a weak 
coupling expansion not supposing a small value for the Higgs coupling, take the 
thermodynamic limit and only then we are allowed to make an expansion at small 
Hq. This procedure can be easily accomplished at least in the main orders of the 
large A expansion if one knows the main state of the gauge system for time-like 
links which has to be perturbated. In the usual ensemble this is a state with all the 
links up or down depending on which n-ality sector we consider in the CE. In the 
F-ensemble there is a competition between one of these states and the state where 
all the time-like spins are flipped inside a volume bounded by S relatively to the 
links outside. The classical action on these two configurations is, correspondingly 



S' = XoNp, + ho{V- 2Q{J:)) , (45) 

Therefore, 



S'' = Ao (iVp„ - 2Nj,) + hoV. (46) 



S^^ - = -2XoN^ + 2hon{E). (47) 

On any finite lattice there always exists such a large Aq that the surface term wins 
and thus the first configuration will represent the main state. Since the PL is not 
flipped in this case in the corresponding triality sector of the CE one finds At = — 1. 
Since, however we have to consider the limit S oo, which is equivalent to f2 — > cxo 
and Nj:, — > oo, we get easily convinced that this state gets metastable when the 
volume term suppresses the surface term in this limit. Hence, it is the second state 
( ^61) which gives the dominating contribution to the thermodynamic limit. Since 
the PL in the origin flips the sign one finds At = 1. 

These qualitative arguments can be made more precise if we consider an effective 
model for the PL which can be obtained in the region Aq >> A„, >> 
Summing over space gauge fields and over Higgs fields we come to the standard 
effective model for the PL which in this case reads (we omit all irrelevant constants) 



^e//= E exp 



(48) 



where x denotes now the sites of the D dimensional lattice. This is the familiar 
Ising model in the external magnetic field with effective couplings defined as 

tanh7 = (tanhAo)^* , tanha = (tanh/io)^*- 

In the usual ensemble ax = a, while in the F-ensemble one has to change sign in 
a if X & n(E). We calculated the correlation function in the F-ensemble and the 
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quantity At in the Ising model both analytically and in MC simulations [O. The 



results found support completely the qualitative arguments we presented above. The 
interpretation of such behaviour we gave in the previous section. 

It is interesting to mention that there is, at least, one obvious possibility which 
leads to a nontrivial value At = —1 even in the thermodynamic limit. This is 
the case of an inhomogeneous magnetic field in the effective Ising model (0). For 
instance, one can take a slowly varying field of the form 

271 , , 

ax = ao cos -j-x, (49) 

where L is the linear size of the system. In such a field one expects that the ground 
state is again a state with spins flipped inside a volume enclosed by an appropriate S 
and, thus A = —1 at large values of 7. At small 7 A = 1 as can be easily shown via 
the high temperature expansion. This implies that in the model with Ising spins in 
the external field (^) there could be a phase transition. With respect to the Higgs 
coupling it means that h{x, t) may be a constant everywhere except one time-slice 
where it has to be of the form (^). Whether such a field could have any physical 
significance in the context of gauge models remains however unclear. 

We turn now to the order parameter A^. Let Ws{C) be a space-like Wilson 
loop in the finite temperature gauge-Higgs model defined in (^) and be a two 
dimensional surface on a dual lattice, fls ^ the corresponding volume. Following the 
above procedure we have 

The temporal part of the Higgs action is not affected by the singular Z{2) gauge 
transformations whereas for the spatial part we get in the F-ensemble 

x,n 

where /i„(x) is defined similarly to (^). Again, we consider the strong and weak 
coupling phases in the gauge coupling A. 

I. As in the strong coupling region A„ << 1. 

Using the general expansion (|40|) one gets up to second order 

<W{C) >= ntanh/i^l + 2L)PctanhAtanh2/i(l -tanh^/i)], (52) 

lec 

where we consider the symmetric case Aq = A„ and ho = hn, for simplicity. Pc is 
the perimeter of the loop C. It gives 

As = l-Oi\'). 
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As in the previous case for At, it is straightforward to see that 

^ = l-A^s.C2(tanhA)^^^ (53) 

which gives Ag — 1 in the — > oo hmit. 

II. As in the weak coupling region A„ >> 1. 

The gauge part gives the following contribution to the Wilson loop for the space 
dimension D = 3 

< W{C) >cx exp [-2Pcie-^^Y] . (54) 

We again analyze the symmetric case as above. The fermionic screening is now 
suppressed as h^"^, i.e. 

< W{C) >oc (tanh/i„)^^. (55) 

The crucial point distinguishing this case from the similar expansion in the electric 

sector is that now we are allowed to expand in small hn since there are no loops 
going around the lattice in space direction which could destroy the convergence. It 
is straightforward to calculate, e.g. < F{T,s) > in leading order of small h 

< F(E,) >= exp [SN^^ + 0{h^)] , (56) 

where 5 ~ 2/i^ tanh A. If we remove the surface back to the pure gauge action, 
then the main contribution in the F-ensemble comes from configurations of the 
gauge fields Sn{x) fiipped in the volume f2(Es) relatively to Sn{x) outside of r2(Ss). 
Since the linking number of the surface E^ and the Wilson loop C is 1, the Wilson 
loop changes sign in the F-ensemble. The corresponding expansions in large A and 
in small hn are converging. Therefore, in the limit of infinite E^ and C all the 
corrections go to zero and we find 

^(E3,C) = -1. (57) 

This indicates perhaps a deconfinement phase transition in the gauge coupling to 
a phase where the kinetic screening dominates the dynamical screening. Above the 
critical point one can detect a unit of Z{2) fiux by the Aharonov-Bohm effect. In 
the main order the critical line is determined as 



exp 



tanh h": 



Our last comment concerns the usage of the CE. One can easily check that Ag 
exhibits essentially the same behaviour in both ensembles, at least in the main 
order of the Higgs coupling which is not affected by the temporal loops. Thus, both 
ensembles should lead to the same phase structure. The use of the CE is more 
relevant, as we believe, for the investigation of the Z{Nc) symmetry in the electric 
sector, in particular to give a proper description of the metastablc states. 



21 



5 Summary 



In this paper we proposed a method which distinguishes between different screening 
mechanisms in gauge theories at finite temperature. This method rehes on the action 
of the order parameter A used earher in the context of zero temperature theory [|1^ 
and is based on the Aharonov-Bohm effect. Various phenomena can be studied 
using this approach. The most notable apphcation is a possibihty for describing the 
deconfinement phase transition in terms of screening mechanisms. We suggested a 
general line for such an application and illustrated it on the example of Z{2) gauge 
spins coupled to Higgs fields. Among other applications it is worth to mention the 
following: 

1) The quantity At can be used to study the problem of domain walls and 
of the stability of the corresponding interfaces in finite temperature theories with 
dynamical matter fields. That this order parameter indeed may show the instability 
of the interface in the thermodynamic limit we showed in Jl^ on the example of the 
3D Ising model. 

2) A very important and long standing problem in the dynamics of gauge systems 
(not only for finite temperature ones) is the difference between the Higgs and the 
fermionic screening. In particular, it is needed to explain the non-confining character 
of the weak interaction. Since the quantity A provides a detailed description of the 
screening in both systems this problem can be formulated rigorously in the present 
approach. 

3) We conjecture here that a similar behaviour of A can also be found in the 
model of Z{Nc) gauge spins coupled to fundamental fermions. The general phase 
structure of the latter model may however differ from the Higgs model studied here. 
Two essential points distinguish the fermion model from the Higgs model. First of 
all, in the limit of large A one can get a model of a free fermion gas which does 
not possess a phase transition. On the other hand, in the limit of strong coupling 
one can find a phase transition in the fermion sector which is related to chiral 
symmetry restoration. It is an interesting opportunity to explore the properties of 
A in both electric and magnetic sectors in order to investigate a phase structure of 
this theory and to get a deeper insight into the connection (if any) between screening 
mechanisms of triality and chiral symmetry restoration. 

4) The application of all previous suggestions and results to SU(Nc) gauge models 
at finite temperature is, in principle straightforward though it is technically much 
more involved. Summarizing we would like to stress that although a final answer to 
the issue of deconfinement in QCD is not available yet, certainly the order parameter 
A should be able to check the scenario proposed here also in the context of full QCD. 

We would like to thank J. Polonyi for interesting discussions during preparation 
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of this work in Vienna and Debrecen. 
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